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Abstract 

It is known how to construct, in a bipartite quantum system, a unique low rank entangled 
| mixed state with positive partial transpose (a PPT state) from an unextendible product basis (a 

UPB), defined as an unextendible set of orthogonal product vectors. We point out that a state con- 
structed in this way belongs to a continuous family of entangled PPT states of the same rank, all 
related by non-singular product transformations, unitary or non-unitary. The characteristic prop- 
erty of a state p in such a family is that its kernel Ker p has a generalized UPB, a basis of product 
vectors, not necessarily orthogonal, with no product vector in Im p, the orthogonal complement 
of Ker p. The generalized UPB in Ker p has the special property that it can be transformed to 
orthogonal form by a product transformation. In the case of a system of dimension 3 x 3, we 
give a complete parametrization of orthogonal UPBs. This is then a parametrization of families 
of rank 4 entangled (and extremal) PPT states, and we present strong numerical evidence that it is 
a complete classification of such states. We speculate that the lowest rank entangled and extremal 
PPT states also in higher dimensions are related to generalized, non-orthogonal UPBs in similar 
ways. 
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^ ! 1 Introduction 

O . 

For a composite quantum system, with two separate parts A and B, the mixed quantum states are 
described by density matrices that can be classified as being either entangled or separable (non- 
entangled). However, there is in general no easy way to classify a given density matrix as being 
separable or not. This problem is referred to as the separability problem, and it has been approached 
in the literature in different ways over the past several years (1). As a part of this discussion there has 
been a focus on a subset of the density matrices which includes, but is generally larger than, the set of 
separable states. This is the set of the so-called PPT states, the density matrices that remain positive 
under a partial matrix transposition, with respect to one of the subsystems, either A or B O. 

Since it is straight forward to establish whether a density matrix is PPT, the separability problem is 
reduced to identifying the subset of entangled PPT states. We refer here to the set of separable states 
as S and the set of PPT states as V, with S C V. These are both convex subsets of the full convex set 
of density matrices, which we denote as V, and in principle the two sets are therefore defined by their 
extremal states. The extremal separable states are the pure product states, and these are also extremal 
states of the set V. Since V is in general larger than <S, it has additional extremal states, and these 
states are not fully known. The problem of finding and classifying these additional extremal states is 
therefore an important part of the problem to identify the PPT states that are entangled. 
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We have in two previous publications studied, in different ways, the problem of finding extremal 
PPT states in systems of low dimensions. In [3] a criterion for extremality was established and a 
method was described to numerically search for extremal PPT states. This method was applied to 
different composite systems, and several types of extremal states were found. In a recent paper [4] 
this study has been followed up by a systematic search for PPT states of different ranks. Series of 
extremal PPT states have there been identified and tabulated for different bipartite systems of low 
dimensions. 

The study in [4] seems to show that the extremal PPT states with lowest rank are somehow special 
compared to the other extremal states. In particular we have found that these density matrices have 
no product vectors in their image, but a finite, complete set of product vectors in their kernel. This 
was found to be a common property of the lowest rank extremal PPT states studied there, for all 
systems with subsystems of dimensions larger than 2. This property relates these states to a particular 
construction, where unextendible product bases, UPBs for short, are used in a method to construct 
entangled PPT states EH0. 

The motivation for the present paper is to follow up this apparent link between the lowest rank 
extremal PPT states and the UPB construction. Our focus is particularly on the rank 4 states of the 
3x3 system. The rank 4 extremal PPT states that we find numerically by the method introduced in 
[4] are related by product transformations to states constructed directly from UPBs. We discuss this 
relation and use it to give a parametrization of the rank 4 extremal PPT states. 

Although a direct application of the (generalized) UPB construction to the lowest rank extremal 
states is restricted to the 3x3 system, the similarity between these states and the lowest rank extremal 
states in higher dimensions indicates that there may exist a generalization of this construction that is 
more generally valid. We include at the end a brief discussion of the higher dimensional cases and 
only suggest that a construction method, and thereby a parametrization, of such states may exist. 

2 An extension of the UPB construction of entangled PPT states 

We consider in the following a bipartite quantum system with a Hilbert space H = Ha ® 'Hb of 
dimension N = NaNb- By definition, a separable state can be written as a density operator of the 
form 



with ph > 0, J2kPk = 1> an d with tpk = 4>k ® Xk as normalized product vectors. The image of 
p, Im p, is spanned by these vectors. The fact that Im p must be spanned by product vectors if p is 
separable is the basis for the UPB construction, which was introduced in Ref. [5], and used there to 
find low-rank entangled PPT states of the 3x3 system. We review here this construction and discuss a 
particular generalization. 

Consider U to be a subspace of % that is spanned by a set of orthonormal product vectors 



which cannot be extended further in % to a set of p + 1 orthogonal product vectors. This defines the 
set as an unextendible product basis (a UPB). Let U be the orthogonal complement to U. The state 
proportional to the orthogonal projection onto U.^, 




(1) 



k 



tpk = 4>k®Xk , k= 1,2,..., p 



(2) 




(3) 
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with a\ = 1/ (N — p) as a normalization factor, is then an entangled PPT state. It is non-separable 
because Im p\ = U 1 - contains no product vector, and it is PPT because pf , the partial transpose of p\ 
with respect to subsystem B, is proportional to a projection of the same form, 

Pi=ai(l-^2^i\, (4) 

with fa = fa ® x*k ■ The vector x% is the complex conjugate of Xk > m the same basis in T~Lb as is 
used for the partial transposition. 

The set of product vectors { fa = fa <g> x% } is a new orthonormal UPB, which generally spans 
a different subspace than the original set { fa = fa (g> Xfc }■ However, it may happen that there 
exists a basis for the Hilbert space T-Lb of the second subsystem in which all the vectors Xk have real 
components. In such a basis the two UPB sets are identical and the state p\ is PPT for the simple 
reason that it is invariant under partial transposition, pf = p\. All the states given as examples in 
Ref. are of this special kind. 

An entangled PPT state p\ defined by this UPB construction is a rather special density operator. 
Ker p\ is spanned by product vectors, while Im p\ contains no product vector. Since p\ is proportional 
to the orthogonal projection onto the subspace Im p\, it is the maximally mixed state on this subspace. 
There is also a symmetry between p\ and pf, such that pf shares with pi all the properties mentioned 
above, and has the same rank N — p, where N = NaNb is the dimension of the Hilbert space and p 
is the number of product vectors in the UPB. 

Implicitly the construction implies limits to the rank of p\. Thus, for a given Hilbert space of 
dimension N = NaNb there is a lower limit to the number of product vectors in a UPB, which 
follows from the requirement that there should exist no product vector in the orthogonal space U ± . 
The corresponding upper bound on the rank to of pi, as discussed in Ref. [4], is given by to < 
N — Na — Nb + 2. There is also a lower bound to > maxjA^, Nb}, which is the general lower 
bound on the rank of entangled PPT states with full local rank [8]. In some special cases there exist 
more restrictive bounds than the ones given here |9). 

For the 3x3 system these two bounds allow only one value m = 4 for the rank of a state p\ 
consructed from a UPB, and for this rank explicit constructions of UPBs exist 0. Also in higher 
dimensions a few examples of UPB constructions have been given |6j. 

The extension of the UPB construction that we shall consider here is based on a certain concept 
of equivalence between density operators previously discussed in iflOl . The equivalence relation is 
defined by transformations between density operators of the form 

92 = a 2 V Pl V^ , (5) 

where a 2 is a positive normalization factor, and V = Va 8) Vb, with Va and Vb as non-singular 
linear operators on Ha and Hb, respectively. The operators p\ and p 2 are equivalent in the sense 
that they have in common several properties related to entanglement. In particular, the form of the 
operator V implies that separability as well as the PPT property are preserved under the transformation 
(f5]). Preservation of separability follows directly from the product form of the transformation, while 
preservation of PPT follows because the partially transposed matrix pf is transformed in a similar 
way as pi, 

pi = a 2 V Pl p V^ , (6) 

with V = Va <8> Vg. If Pi and pf are both positive then the transformation equations show explicitly 
that the same is true for p 2 and p 2 p . Furthermore, since the operators V and V are non-singular, the 
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ranks of p\ and p 2 are the same, and also the ranks of pf and p 2 p . The same is true for the local 
ranks of the operators, which are the ranks of the reduced density operators, defined with respect to 
the subsystems A and B. Finally, if p\ is an extremal PPT state, so is p 2 . 

Let us again assume p\ to be given by the expression ©. Since the product operator V is an 
invertible mapping from Im p\ to Im p 2 , and since Im p\ contains no product vector, there is also no 
product vector in Im p 2 , hence p 2 is entangled. Similarly, the product operator {V^)^ 1 is an invertible 
mapping from Ker p\ to Ker p 2 , and it maps the UPB in Ker p\, eq. (121), into a set of product vectors 
in Ker p 2 , 

^ = ((yjrVfc) ® ((^Jrtfc) , k=i,2,..., P . (7) 

If the operators Va and Vb are both unitary, then this is another UPB of orthonormal product vectors, 
and p 2 is proportional to a projection, just like p\. More generally, however, we may allow Va and 
Vb to be non-unitary. Then the product vectors ip' k in Ker p 2 will no longer be orthogonal, but p 2 is 
nevertheless an entangled PPT state. It has the same rank as p\ , but is not proportional to a projection. 

Since the normalization of the density operators p\ and p 2 is taken care of by the normalization 
factors a\ and a 2 , we may impose the normalization condition det Va = det Vb = 1, which defines 
the operators as belonging to the Special Linear (SL) groups on %a and %b- We will say then that the 
two density operators p\ and p 2 , related by a transformation of the form (f5]), are SL ® SL equivalent, 
or simply SL equivalent. 

The above construction motivates a generalization of the concept of a UPB, where we no longer 
require the product vectors to be orthogonal. This generalization has also previously been proposed 
in the literature [7]. In the following we will refer to an unextendible product basis of orthogonal 
vectors as an orthogonal UPB. A more general UPB is then a set of product vectors that need not 
be orthogonal (need not even be linearly independent), but satisfies still the condition that no product 
vector exists in the subspace orthogonal to the set. The UPB defined by © is a special type of 
generalized UPB, since it is SL equivalent to an orthogonal UPB. More general types of UPBs exist, 
and they are in fact easy to generate, since an arbitrarily chosen set of k product vectors is typically a 
generalized UPB, in the above sense, when k is sufficiently large. However, if it is not SL equivalent 
to an orthogonal UPB, then we have no guarantee that there will be any entangled PPT state in the 
subspace U orthogonal to the generalized UPB. 

3 Parametrizing the UPBs of the 3x3 system 

We focus now on the orthogonal UPBs in the 3x3 system, which must have precisely 5 members. In 
fact, for any given set of 4 product vectors <pk (g> Xk, there exists a product vector (j)®x orthogonal to 
all of them, for example with cj)\ _L _L <p 2 and %3 -L X -L X4- And with 6 members in an orthogonal 
UPB, it would define a rank 3 entangled PPT state, which is known not to exist (H. 

The general condition for 5 product vectors to form an orthogonal UPB in the 3x3 system was dis- 
cussed in Ref. 0. The condition implies that for any choice of three product vectors from the set, the 
first factors 4>h are linearly independent and so are the second factors \k- The orthogonality condition 
further implies that if the product vectors are suitably ordered, there is a cyclic set of orthogonality 
relations between the factors of the products, of the form 

01 ± <p 2 _L 4> 3 _L 04 _L <p 5 _L fa , xi -L X3 -L X5 -L X2 -L Xa -L Xi • (8) 

In Fig. [TJ the situation is illustrated by a diagram composed of a pentagon and pentagram, where each 
corner represents a product vector. Each pair of vectors is interconnected by a line showing their 
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orthogonality. A solid (blue) line indicates orthogonality between (ft states (of subsystem A) and a 
dashed (red) line indicates orthogonality between x states. As shown in the diagram, precisely two A 
lines and two B lines connect any given corner with the other corners of the diagram. 




Figure 1: Diagrammatic representation of the orthogonality relations in a 5-dimensional UPB of the 3x3 sys- 
tem. The corners of the diagram represent the product vectors of the UPB, and the lines represent orthogonality 
between pairs of states. There are two types of orthogonality, represented by the solid (blue) lines and the 
dashed (red) lines. The solid lines represent orthogonality between the vectors of the products that belong to 
subsystem A and the dashed lines represent orthogonality between the vectors belonging to subsystem B. 

Introducing a complete set of orthonormal basis vectors ccj in Ha> we write 

3 

<ftk = Ujk <Xj i k = 1,2,3,4, 5 . (9) 

3=1 

We may choose, for example, a% proportional to cft\ and «2 proportional to <ft>2- If we multiply each 
basis vector aj by a phase factor uij, and each vector 4>k by a normalization factor N^, we change 
the 3x5 matrix Ujk into uj^ l NkUjk- It is always possible to choose these factors so as to obtain a 
standard form 

/10a b 0\ 
u= 1 1 a , (10) 
\0 b -a 1/ 

with a and b as real and strictly positive parameters, and with the vectors (ftk not normalized to length 
1. A similar parametrization of the vectors of subsystem B with orthonormal basis vectors gives 

3 

Xk = ^ l v jk P j , k = 1,2,3,4,5, (11) 

3=1 

and a standard form 

/l d c\ 
v = 1 1 c , (12) 
\0 -c 1 d) 

with two more positive parameters c and d. Thus, an arbitrary orthogonal UPB is defined, up to unitary 
transformations in Ha and 7-(.B,by four continuous, positive parameters a, 6, c, d. 
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Note that, for a given UPB, the parameter values are not uniquely determined, since the above 
prescription does not specify a unique ordering of the 5 product vectors within the set. Any permuta- 
tion that preserves the orthogonality relations pictured in Fig. [T] will generate a new set of values of 
the parameters that define the same UPB. These permutations form a discrete group with 10 elements, 
generated by the cyclic shift k — >• k + 1, and the reflection k — > 6 — k. 

Given the orthonormal basis vectors ay in Ha and /3j in Hb, we may think of the four positive 
parameters a, 6, c, d as defining not only one single orthogonal UPB, but a continuous family of gener- 
alized UPBs that are SL equivalent to this particular orthogonal UPB. The parameter values defining 
one such family may be computed from any UPB in the family via SL invariant quantities, in the 
following way. Given the product vectors <pk (g> Xk for k = 1, 2, 3, 4, 5, not necessarily orthogonal, we 
introduce expansion coefficients as in (© and arrange them as column vectors 



Uk = \ u 2 k • (13) 




Then we introduce the following quantities, 

det(uiU2'U4) det(ui«3ti5) 2 

1 det(niM2^5) det(ni^3U4) 

^ _ det(uiU2U 3 ) det(u 2 U4«5) _ b 2 

det(uiU2U4 : ) det(u2U3Us) a 2 ' 

where the values to the right are determined from the parametrization (TfOl of the orthogonal UPB 
defining the family. Similarly, we define 



det(viV2^3) detfviv^v^) 

S3 = = c 

det(uiU2U5) det(i> 1V3V4) 



2 



s4 _ det^ii^) det(f 2 ^3^4) _ (P_ ^ 
det(^iV2^3) detfazViVs) c 2 

The quantities si, S2, S3, S4 defined in terms of 3 x 3 determinants are useful because they are invariant 
under SL transformations as in (f7]), and in addition they are independent of the normalization of the 
column vectors and v^. Obviously, many more similar invariants may be defined from 5 product 
vectors, but these four invariants are sufficient to characterize a family of UPBs that are SL equivalent 
to an orthogonal UPB. 

There exists a less obvious further extension of the set of invariants. In fact, there are always 6 
vectors that can be used to define invariants, since in addition to the 5 linearly independent product 
vectors of the UPB, the space spanned by these will always contain a 6th product vector. In the case 
of an orthogonal UPB, given by the parameters a, b, c, d, we have found (by means of a computer 
algebra program) explicit polynomial expressions for the components of the one extra product vector. 
We have checked, both analytically and numerically, that the existence of exactly 6 product vectors 
is a generic property of a 5 dimensional subspace of the 3x3 dimensional Hilbert space %. This 
number of product vectors is also consistent with the formula discussed in H, which specifies more 
generally, as a function of the dimensions, the number of product vectors in a subspace of %. For 
an orthogonal UPB in the 3x3 system the 6th vector is singled out because it is not orthogonal to the 
other vectors, but for a non-orthogonal UPB there is no intrinsic difference between the 6 vectors of 
the set, which should therefore be treated on an equal footing. 
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For a UPB that is SL equivalent to an orthogonal UPB there are strong restrictions on the values of 
invariants of the above kind, since they are all rational functions of the four real parameters a, b, c, d. 
In particular, they must all take real values. A given choice of four invariants, as in (fT4l and (fl~5T ), 
is sufficient to define the parameter space for the equivalence classes of these UPBs. But since the 6 
product vectors listed in any order define the same UPB, and the same PPT state, there is a discrete 
set of 6! = 720 symmetry transformations that introduce identifications between points in the corre- 
sponding parameter space. As we shall see below, the requirement that all four invariants si, S2, S3, S4 
should be positive allows 60 different orderings from the total of 720. 

One should note that for a generalized UPB consisting of 5 randomly chosen product vectors the 
invariants will in general be complex rather than real, and it is not a priori clear that four invariants 
are sufficient to parametrize the equivalence classes of random UPBs. 

4 Classifying the rank 4 entangled PPT states 

We have in [4 ] described a method to generate PPT states p for given ranks (m, n) in low-dimensional 
systems, with m = rankp and n = rankp p . By repeatedly using this method with different initial 
data we have generated a large number of different PPT states of rank (4, 4) in the 3x3 system. They 
are all entangled PPT states, and as a consequence they are extremal PPT states. This follows from 
the fact that if they were not extremal they would have to be convex combinations involving entangled 
PPT states of even lower ranks, and such states do not exist. 

The remarkable fact is that every one of these states has a UPB in its kernel which is SL equivalent 
to an orthogonal UPB, and the state itself is SL equivalent to the state constructed from the orthogonal 
UPB. We regard our numerical results as strong evidence for our belief that the four real parameters 
which parametrize the orthogonal UPBs give a complete parametrization of the rank 4 entangled PPT 
states of the 3x3 system, up to the SL (or more precisely SL cg> SL) equivalence. We will describe 
here in more detail the numerical methods and results that lead us to this conclusion. 

Assume p to be an entangled PPT state of rank (4, 4), found by the method described in [4]. The 
question to examine is whether it is SL equivalent to an entangled PPT state defined by the orthogonal 
UPB construction. We therefore make the ansatz that it can be written as p = p2 = a^V p\V\ where 
pi is defined by a so far unknown orthogonal UPB, parametrized as in (TTOb and (fT2l . and where 
the transformation V is of product form, V = Va <8>Vb- We consider how to compute the product 
transformation V, assuming that it exists. The fact that we are able to find such a transformation for 
every (4, 4) state is a highly non-trivial result. 

Given p, the first step is to find all the product vectors in Ker p. We solve this as a minimization 
problem: a normalized product vector ip = <p <8> X with ptp = is a minimum point of the expectation 
value tp'pip. Details of the method we use are given in Ref. [4]. Empirically, we always find exactly 
6 such product vectors ipf, = <pk <8> k = 1, 2, . . . , 6, any 5 of which are linearly independent and 
form a UPB, typically non-orthogonal. 

Although the numbering of the 6 product vectors is arbitrary at this stage, we compute the in- 
variants si, S2, S3, S4, substituting for uu and Xk for u&, with k = 1, 2, . . . , 5. As shown by the 
previous discussion all the four invariants have to be real, otherwise no solution can exist. A random 
UPB has complex invariants, and the empirical fact the invariants are always real for a UPB in Ker p 
where p is a rank (4, 4) entangled PPT state, is a non-trivial test of the hypothesis that such a UPB can 
be transformed into orthogonal form. 

It is not sufficient that the invariants are real. As shown by the expressions (fT4l and (fl3T ) there 
has to exist an ordering of the product vectors where all the four invariants take positive values. The 
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signs of the invariants will depend on the ordering of the product vectors, and most orderings produce 
both positive and negative invariants. For the rank (4, 4) density matrices that we have constructed, it 
turns out that it is always possible to renumber the 5 first vectors in the set in such a way that all four 
invariants become positive. This is a further non-trivial test of our hypothesis. 

There are in fact, in all the cases we have studied, precisely 10 of the 5! permutations of the 5 
vectors that give positive values of the four invariants. This means that such an ordering is unique 
up to the symmetries noticed for the diagram in Fig. [T] However, there is a further symmetry, since 
the reordering which gives positive invariants works for any choice of the 6th vector of the set. The 
possible reorderings of all 6 product vectors which preserve the positivity of the invariants therefore 
define a discrete symmetry group with altogether 6 x 10 = 60 elements, which defines mappings 
between different, but equivalent, representations of the UPB in terms of the set of four real and 
positive invariants. The corresponding parameter transformations are given in the Appendix. 

Assume now, for a given rank (4, 4) state, that a "good" numbering has been chosen for the 6 
product vectors tp k = 4> k ® Xk in the corresponding UPB, so that the four invariants defined by the 
first 5 vectors are all real and positive. The problem to be solved is then to find the transformation 
that brings the UPB into orthogonal form. This means to find 3x3 matrices C and D such that 
4>k = N' k Cuk and Xk = N' k Dv k for k = 1, 2, . . . , 5, with unspecified normalization constants N' k 
and N' k '. Here the vectors u k and v k belong to an orthogonal UPB as given by the equations (fTOl) 
and (fl~2~l) . and these vectors are all known at this stage, because the invariants s%, S2, S3, S4 determine 
the parameters a, b, c, d. The transformation matrices C and D correspond to and Vg in ©. The 
condition for two vectors (f) k and Cu k to be proportional is that their antisymmetric tensor product 
vanishes, hence we write the following homogeneous linear equations for the matrix C, 

4>k A (Cu k ) = 4> k ® (Cu k ) - (Cu k ) ® 4> k = , k = 1,2,..., 5 . (16) 

Since the antisymmetric tensor product <j) k A (Cu k ) has, for given k, 3 independent components, there 
are altogether 15 linear equations for the 9 unknown matrix elements Cy. We may rearrange the 3x3 
matrix C as a 9 x 1 matrix C and write a matrix equation 

MC = , (17) 

where M is a 15 x 9 matrix. This equation implies that (M'M)C = 0. The other way around, the 
equation (M^M)C = implies that (MC)t(MC) = C^(M^M)C = and hence MC = 0. Thus 
we may compute the matrix C as an eigenvector with zero eigenvalue of the Hermitean 9x9 matrix 
M'M. The matrix D is computed in a similar way. 

It is a final non-trivial empirical fact for the (4, 4) states we have found, that solutions always exist 
for the matrices C and D, whenever the ordering of the 6 product vectors tp k = <p k <g> \k is such that 
the invariants s\, S2, S3, S4 are positive. 

The result is that every rank (4, 4) state of the 3x3 system which we have found in numerical 
searches [4] can be transformed into a projection operator with an orthogonal UPB in its kernel. We 
have also checked the published examples of entangled PPT states of rank (4,4), which are based 
on special constructions E |U [TT] [12]], and have got the same result for these states. The explicit 
transformations have been found numerically by the method discussed above, and in all cases the four 
parameters a, b, c, d have been determined, with values that are unique up to arbitrary permutations of 
product vectors from the 60 element symmetry group. 
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5 Summary and outlook 



The main result of this paper is a classification of the rank 4 entangled PPT states of the 3x3 system. 
We find empirically that every state of this kind is equivalent, by a product transformation of the form 
SL ® SL, to a state constructed from an orthogonal unextendible product basis. We refer to this type 
of equivalence as SL equivalence. We have shown how to parametrize the orthogonal UPBs by four 
real and positive parameters, and we have described how permutations of the vectors in the UPB give 
rise to identifications in the four-parameter space. 

The concept of SL equivalence of states and of product vectors leads to a generalization of the 
concept of unextendible product bases so as to include sets of non-orthogonal product vectors, and 
further to the concept of equivalence classes of generalized UPBs that are SL equivalent to orthogonal 
UPBs. Thus, the parametrization of the orthogonal UPBs by four positive parameters is at the same 
time a parametrization of the corresponding equivalence classes of generalized UPBs. 

We have described a method for checking whether a given rank 4 entangled PPT state in the 3x3 
system is equivalent, by a product transformation, to a state constructed from an orthogonal UPB. It 
is a highly non-trivial result that all the rank-four entangled states that we have produced numerically, 
and all states of this kind that we have found in the literature, are SL equivalent to states that are 
generated from orthogonal UPBs. This we take as a strong indication that the parametrization of the 
UPBs in fact gives also a parametrization of all the equivalence classes of rank 4 entangled PPT states 
of the 3x3 system. 

Apart from the pure product states, the rank 4 entangled PPT states are the lowest rank extremal 
PPT states among the 3x3 states that we have found in numerical searches, as reported on in [4J. The 
property of such a state, that it has a non-orthogonal UPB in its kernel, which means that there is a 
complete set of product vectors in Ker p and no product vector in Im p, is shared with the lowest rank 
extremal PPT states of the other systems that we have studied, of dimensions different from 3x3. This 
has led us to conjecture that this is a general feature of the lowest rank extremal PPT states, valid 
also in higher dimensional systems [4], and to speculate that there may exist a generalization of the 
construction used for the 3x3 system in terms of orthogonal UPBs and SL transformations, which can 
be applied in the higher dimensional systems. 

In higher dimensions the orthogonality condition is harder to satisfy, and therefore another condi- 
tion may take its place as the defining characteristic of a special subset of extremal states from each 
SL equivalence class. This hypothetical new condition may involve the full set of product vectors in 
the subspace, rather than an arbitrarily selected subset as in the definition of the orthogonal UPBs. 
To examine this possibility, with the aim of parametrizing the lowest rank extremal PPT states more 
generally, we consider an interesting task for further research, and we are currently looking into the 
problem. 
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A Equivalent orderings of the 6 product vectors 



Assume that the sequence of product vectors ipk = 4>k ® Xfc, k = 1, 2, 3, 4, 5, in this order, is charc- 
terized by parameter values a, b, c, d, as computed from the invariants s\, S2, S3, S4. It is convenient 
here to replace the parameters a, b, c, d by a = a 2 , (3 = b 2 , 7 = c 2 , S = d 2 . 

Then the cyclic permutation tp^ h-> with tpi = tp 5 and ^ = ipk-i f° r & = 2, 3,4, 5 corre- 
sponds to the following parameter transformation, which is periodic with period 5, 



a 



1 + Q ' 

P 



a(l + a) 
1 

7(1 + 7 + 6) 
(5(7 + S) 



7 = 37T , (18) 

7 + 

5 = 



The inversion ipi i-> ^1 = ipi, ipk •->■ ^fc = V^-fc f° r & = 2, 3, 4, 5 corresponds to the parameter 
transformation a = a, 7 = 7, 



~ = a(l + a) 



« = ^7-^ . (19) 


Let ^6 = 4>& ® Xe t> e tne 6th product vector in the 5 dimensional subspace spanned by the above 
5 product vectors. Then the sequence tpi = tpQ, tp2 = ^3 = ^3> = ^5 = ^2 corresponds 
to the parameter transformation a = 7, 7 = a, 

= (3(1 + 7)((a + + 5) + 6) 



a(l + a + /3)<5 + (1 + a)(a + /3)(1 + 7) ' 

c = (l + a)(^+(a + /3)7(l+7 + ^)) f2m 
(l + a + (l + a + /?)( 7 + <5))<5 ' 1 ; 

It is not easy to see by looking at the formulae that this parameter transformation is its own inverse. 

Altogether, these transformations generate a transformation group of order 60 (with 60 elements), 
isomorphic to the symmetry group of a regular icosahedron with opposite corners identified. Equiv- 
alently, it is the group of proper rotations of the icosahedron, excluding reflections. The icosahedron 
has 12 corners, and we may associate the 6 product vectors with the 6 pairs of opposite corners. 
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